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Nous nous interessons a la reponse mecanique au cisaillement d’un milieu granulaire modele 
bidimensionel au-dela de la transition de blocage appelee aussi ’’jamming”. Tout d’abord, nous 
developpons le dispositif experimental et nous combinons des techniques de suivi de particules et 
de photoelasticte afin de mesurer l’etat de deformation et l’etat de contrainte a l’echelle du grain. 
Ensuite, nous mettons en place un intrus capable d’extension radiale (un ’’ballon” 2D) afin de pou- 
voir cisailler l’empilement granulaire tout en conservant une geometrie axisymetrique. Nous sondons 
l’apparition des contraintes a l’echelle du grain pour des amplitudes de deformation inferieures a 
10 -2 et pour une gamme de fraction volumique evoluant de 2% de part et d’autre de la transition 
de blocage. Nous montrons ainsi que cette reponse mecanique induit des contraintes de cisaillement 
mais aussi des contraintes normales. De plus, nous identifions un regime elastique ou les contraintes 
normales et les contraintes de cisaillement evoluent non-lineairement avec la deformation de cisaille¬ 
ment. Enfin, nous explicitons la relation entre l’apparition de cette non-linearite et la transition de 
blocage et nous determinons les relations constitutives de l’empilement. 

We investigate experimentally the mechanical response to shear of a 2D packing of grains across 
the jamming transition. First, we develop a dedicated experimental setup, together with tracking and 
photoelastic techniques in order to prepare the packing in a controlled fashion and to quantify the 
stress and strain tensors at the grain scale. Second, we install a inflating probe (a 2D ’’balloon”), 
which shears the packing with a cylindrical symmetry. We probe experimentally stresses and strains 
for strain amplitudes as low as 10 -3 and for a range of packing fractions within 2% variation around 
the jamming transition. We observe not only that shear strain induces shear stresses, but also normal 
stresses. Moreover, we show that both shear and normal stresses behave nonlinearly with the shear 
strain. Finally, we show by scaling analysis that the constitutive laws are controlled by the Jamming 
transition. 

PACS numbers: 45.70.-n 83.80.Fg 


I. INTRODUCTION 

Comprendre les proprietes mecaniques des empile¬ 
ments denses de particules athermiques telles que les 
grains, les mousses ou bien encore les emulsions demeure 
un defi aussi conceptuel que pratique. Les fluctuations, 
le desordre et l’anisotropie controlent la mecanique de 
ces systemes et contrecarrent la plupart des tentatives 
de determination des lois constitutives de ces materiaux. 
Un progres considerable a ete effectue lors de Intro¬ 
duction d’un modele simple de spheres molles non frot- 
tantes, contenant ainsi un degre minimal de complica¬ 
tion. Dans ce modele, un empilement rigide de particules 
compressees les unes contre les autres perd sa stabilite 
mecanique lorsque la fraction d’empilement <f> decroit, en 
dessous d’une fraction particuliere <f> = </>j, qui depend 
du protocole de preparation de cet empilement et de la 
taille du systeme fT |. A ce point particulier, la pression de 
confinement tend vers 0 et les deformations des particules 
disparaissent [2H5]. Les modules elastiques se comportent 
en lois d’echelles avec la distance a ce point. Lorsque 
Ton approche cette transition, le materiau devient alors 
fragile fir et sa reponse lineaire est dominee par les fluc¬ 
tuations a basses energies [7j. Des travaux anterieurs son- 
dant les proprietes mecaniques des empilements de grains 


mettaient en avant que la friction [8| et les chaines de 
forces El HQ] sont a l’origine la nature solide de l’empi- 
lement. Recemment, de nombreux efforts theoriques ont 
montre que la nature marginale de la transition de blo¬ 
cage est l’ingredient essentiel a l’origine des comporte- 
ments non triviaux de la structure de l’empilement et de 
ses proprietes mecaniques et dynamiques [TTJ 12] . II s’agit 
de savoir si les memes comportements apparaissent dans 
les systemes reels. Plusieurs campagnes experimentales 
menees au-dela de la transition de blocage, ont ca- 
racterise les proprietes structurelles et dynamiques des 
systemes granulaires Q3HIE des mousses m et des 
emulsions Q3- II existe aussi quelques etudes qui ont 
determine la rheologie sous la transition de blocage dans 
des experiences de grains vibres [HUH], de mousses ESI, 
et d’emulsions mum. Cependant le lien direct avec la 
transition de blocage n’est pas clairement etabli. En par¬ 
ticulier, la pertinence de la reponse lineaire proche de 
la transition reste au centre des debats j23l - f25| . Pour 
une deformation de cisaillement finie 7, les effets non- 
lineaires deviennent predominants [261428] et la reponse 
mecanique du systeme n’est plus uniquement decrite par 
la geometrie de l’empilement. Enfin, dans le cas des 
experiences avec des grains, les etudes numeriques et 
theoriques de spheres molles ignorent systematiquement 
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les effets de dilatance, c’est a dire l’accroissement de vo¬ 
lume ou de pression sous une deformation de cisaille- 
ment p9ff3T] . 

Dans cet article, qui presente le detail d’une lettre pu- 
bliee recemment |32j, n °us presentons la premiere me- 
sure experimentale de la reponse elastique d’un empi- 
lement 2D de grains a travers la transition de blocage. 
Nous appliquons une deformation de cisaillement inho¬ 
mogene en gonflant un intrus au centre d’une couche bi- 
disperse de grains frottants. Nous determinons le champ 
de deplacement des grains et le reseau de forces de 
contact a partir de mesures photoelastiques et des tech¬ 
niques de suivi de particules, et nous calculons le tenseur 
des deformations et le tenseur des contraintes a l’echelle 
du grain. Les relations constitutives obtenues a partir 
des courbes parametriques des invariants de ces deux 
tenseurs montrent que l’elasticite lineaire ne s’applique 
pas. La dilatance est essentielle et, au-dela du blocage, 
une attenuation de la variation de contrainte de cisaille¬ 
ment se produit pour de faibles deformations. L’elasticite 
lineaire est retrouvee pour de grandes deformations, a 
partir d’une deformation critique q c , qui evolue avec la 
distance au blocage et qui disparait a f>j. Nous recalcu- 
lons enfin les profils de deformations issus des relations 
constitutives et nous montrons qu’ils correspondent aux 
profils experiment aux. 


II. DIPOSITIF EXPERIMENTAL ET 
PROTOCOLE 

Le dispositif est adapte de celui utilise par m ns. 
Une couche bi-disperse de 8166 disques photoelastiques 
de diametre 4 et 5 mm est confinee dans une cellule 
rectangulaire. Un des murs est un piston mobile qui 
permet de controler precisement la fraction surfacique 
(j). Les grains reposent sur une plaque en verre qui peut 
etre vibree a une amplitude de 1 cm a une frequence 
de 10 Hz perpendiculairement a la direction du piston. 
L’intrus gonflable est une entretoise en bronze, equipee 
de 9 pistons radiaux, entouree d’un joint torique de 
diametre 2 77 = 26.3 mm et connectee au reseau d’air 
comprime. Lorsque la pression d’air augmente dans 
1’intrus, les pistons radiaux viennent appuyer sur le joint 
torique, assurant une deformation radiale uniforme, 
jusqu’a 2(77 + a) = 28.5 mm. Quand l’alimentation 
en air est coupee, l’elasticite du joint torique permet 
a celui-ci de reprendre sa forme initiale. Le taux de 
dilatation de l’intrus est a* = a /77 £ [1 — 10]%. La 
specificite de ce chargement localise reside dans le fait 
qu’il sollicite le milieu granulaire selon une compression 
radiale par rapport au centre de l’intrus mais aussi selon 
une extension orthoradiale. 

En variant l’amplitude de la sollicitation et la fraction 
volumique d’empilement, nous enregistrons la reponse 
mecanique selon un protocole precis. Tout d’abord, 
nous introduisons l’intrus au centre de l’empilement a 
la fraction d’empilement la plus basse. Ensuite, nous 


comprimons l’empilement a l’aide du piston mural 
jusqu’a un etat fortement bloque de l’empilement tout 
en vibrant la plaque inferieure (se referer a [151 pour 
des details supplementaires). Nous arretons alors la 
vibration et commencons l’acquisition des images tout 
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Figure 1. Dispositif Experimental, (a) 8166 disques 
photoelastiques sont confines dans un cadre rectangulaire, 
dont la surface est controlee par un piston lateral. Au centre, 
un ’’ballon” 2D est gonfle par un dispositif pneumatique, 
ce qui permet de creer un cisaillement inhomogene au sein 
de l’assemblee. La plaque inferieure (gris clair) sur laquelle 
reposent les grains peut vibrer horizontalement et est uti- 
lisee pour homogeneiser l’assemblee entre les differents tests, 
(b) Test de compression uni-axiale sur un seul grain, les 
donnees brutes (ronds bleus) correspondent a la force F (N) 
en fonction du deplacement S (mm) et un ajustement (ligne 
rouge pleine) de la forme F — A6 b , avec A = 13.0 ± 1.0 et 
b — 1.8=b0.1. (c) Photo de (gauche) donnees brutes entre pola- 
riseurs croises et (droite) d’une reconstruction de l’assemblee 
de grains et du reseau de forces. La taille et la couleur des liens 
correspondent a l’intensite de la force normale entre chaque 
grain. 

Experimental Setup, (a) 8166 photoelastic disks are confi¬ 
ned in a rectangular frame, whose surface is tuned by a side 
piston. At the center, a 2D ’’balloon” is inflated pneumati¬ 
cally, thus inducing an inhomogeneous shear in the packing. 
The bottom plate (light gray), on which rest the grains, can 
be vibrated horizontally and is used to homogenize the packing 
between the different mechanical tests, (b) Uniaxial compres¬ 
sion test on a single grain. The raw data (blue disks) indicates 
the force-displacement curve. The red solid line corresponds to 
a fit to F = AS b , with A = 13.0 =b 1.0 and b = 1.8 d= 0.1. (c) 
Picture of (left) the raw cross-polarized pictures and of (right) 
the reconstructed grain packing and force network. The width 
and color of the links correspond to the intensity of the normal 
force between two grains. 
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en augmentant la taille de l’intrus par pas de 1.5% 
environ. A la fin, nous laissons l’intrus reprendre sa taille 
initiale, reenclenchons la vibration, puis decomprimons 
legerement l’empilement et recommencons le cycle de 
mesure. La vibration permet de re-homogeneiser l’etat de 
contrainte dans rempilement pendant les changements 
de fraction volumique, tout en maintenant une structure 
d’empilement identique El US]. 

Le processus de calibration des grains photoelastiques 
se fait independamment via un essai de compression 
uni-axiale en utilisant une machine de traction uni-axiale 
(Instron 5965) equipee d’une cellule d’effort de 100 N. 
Nous observons que la courbe d’essai est non-lineaire, 
avec un exposant b = 1.8 =b 0.1 (Fig. [TJd) , ce qui est 
coherent avec des etudes precedentes [33J T En revanche, 
ce resultat est en disaccord avec la theorie du contact 
ideal [34], qui predit une loi de contact lineaire avec une 
raideur de k = 1 N/mm pour notre systeme. Dans la 
suite, on normalisera les tenseur de contraintes 2D par 
cette raideur. 

Les grains photoelastiques sont eclaires par le dessous 
de la cellule avec une large source de lumiere polarisee 
et uniforme (PHLOX). Les images sont enregistrees a 
l’aide d’une camera CCD haute resolution (2048 x 2048 
pixels), conduisant a une resolution spatiale de 100/im. 
L’acquisition de l’information des positions des grains 
et du signal photoelastique est realisee a l’aide d’une 
roue, equipee alternativement d’un polariseur croise, 
montee sur un moteur pas-a-pas (se referer a na Eg 
pour plus de details). De ces images, nous extrayons les 
positions des grains, que nous pouvons suivre pour ob- 
tenir leurs trajectoires. Ensuite, nous utilisons le champ 
de positions des grains pour construire la triangulation 
de Delaunay et la tesselation de Voronoi que nous 
combinons avec le signal photoelastique (Fig. [IJ} gauche) 
afin d’estimer les forces normales et tangentielles entre 
les grains (Fig. [TJ^ droite). Nous avons alors acces a un 
certain nombre d’observables a l’echelle du grain, telles 
que le nombre de contacts qu’un grain possede avec ses 
voisins ainsi que les tenseurs des deformations et des 
contraintes. Dans la suite, nous nous concentrons sur 
le nombre de contact, l’etat de contrainte et l’etat de 
deformation qui sont des parametres pertinents pour 
obtenir le comportement macroscopique de l’empilement. 


III. ETAT DE CONTRAINTE ET ETAT DE 
DEFORMATION A L’ECHELLE DU GRAIN 

A partir du champ de deplacement U des grains, nous 
calculons le tenseur des deformations e a l’echelle du 
grain [361438] . Pour cela, nous evaluons le tenseur des 
deformations a partir du champ de deplacements des 
grains voisins. Cette formulation permet d’eviter d’ef- 
fectuer un calcul par differences finies [38] . Les compo- 
santes du tenseur de deformation sont estimees a partir 
du champ de deplacement U du grain i et des normales 


aux cotes de sa cellule de Voronoi 5, definies par sa posi¬ 
tion par rapport a ses voisins (Fig.[2^). Ainsi, le gradient 
moyen du champ de deplacement de U correspondant au 
grain i dans S peut etre defini par : 

(VC/) = 1 JJvudS. (1) 

En utilisant le theoreme de la divergence, on obtient 

(VC/) = ^ju- kdL, (2) 

oil L est 1’arete correspondant a la cellule de Voronoi et k 
est la normale sortante de cette arete (Fig. iv En faisant 
l’hypothese que la cellule de Voronoi est un polygone, on 
trouve 

( VC/ ) = |X C/ - fe ’ ( 3 ) 

U L 

ou R est le rayon du grain i considere. Le tenseur des 
deformations e correspond a la partie symetrique de ce 
gradient du champ de deplacement : 

e = \ ((VC/) + (VC/)t) . (4) 

Le principal atout de cette methode est de ne pas utiliser 
de differences finies, ce qui reduit considerablement le 
bruit du calcul. 

Pour calculer le tenseur des contraintes cr a l’echelle 
du grain i, nous combinons la tesselation et les forces 


Figure 2. Determination des tenseurs de contrainte et 
de deformation, (a) Le grain central a un deplacement U et 
les cotes de sa cellule de Voronoi (en rouge) sont definis par 
les vecteurs ki. La somme du produit tensoriel de U avec ki 
donne ainsi le gradient de deplacement, (b) Pour le tenseur 
des contraintes, les vecteurs ki sont combines avec les forces 
de contact fi (en vert). 

Measurement of the stress and strain tensors, (a) The 

inner grain has a displacement U and the edges of its Vo¬ 
ronoi cell (in red) are defined by the vectors ki. The sum of 
the tensorial product between U and ki is the displacement 
gradient, (b) As for the stress tensor, the vectors ki are com¬ 
bined with the contact forces fi (in green). 
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d’interactions du grains i avec ses voisins selon EH [36] : 

<T “|XA‘A’ (5) 

ou fjj est la force de contact entre le grain i et un de ses 
voisins j tandis que est le rayon vecteur oriente du 
centre du grain i vers le centre du grain j (Fig. [2^). 

Apres avoir verifie que ces tenseurs part agent les 

memes vecteurs propres [39], nous restreignons l’analyse 

au premier et second invariants de chaque tenseur. Ainsi, 
la dilatation est definie par 

e = 2 (6) 

k 

la pression est definie par 

p = -\^2^kk, (7) 

k 

la deformation de cisaillement est definie par 



et la contrainte de cisaillement est definie par 



oil Sij represente le symbole de Kronecker. Dans la suite, 
P et r sont normalises par la raideur du contact k = 
1 N/mm et 1’unite de longueur est le diametre des petits 
grains 8 = 4 mm. Le tenseur des contraintes et le tenseur 
des deformations sont respectivement mesures avec une 
resolution de 10 -4 et 10 -3 . 

IV. ETAT DE L’EMPILEMENT NON 
SOLLICITE 

Pour chaque fraction d’empilement, avant le gonfle- 
ment de l’intrus, le systeme est caracterise par un etat 
initial, avec des chaines de force se repartissant dans 
l’ensemble du materiau. La particularity de l’empile- 
ment bloque reside dans le fait que toutes les forces 
sont repulsives, la rigidite de l’assemblee etant donnee 
par une pression de confinement. Un tel etat bloque, 
qui a auparavant ete etudie en detail m , est statisti- 
quement homogene. Le nombre de contact moyen zq est 
essentiellement constant pour les faibles fractions sur- 
faciques (Fig. [3^i) . Aux fractions surfaciques moyennes, 
il presente une singularity a partir de laquelle il aug- 
mente non-lineair ement. On identifie cette singularity 
avec la transition de blocage a la fraction surfacique 
4>j = 0.8251 ± 0.0009. Il n’est pas surprenant d’obser- 
ver une valeur non nulle Zq sous <f>j : quand la vibra¬ 
tion est arretee, la structure est figee brutalement et des 


forces residuelles demeurent a cause de la friction entre les 
grains et la paroi inferieure. La croissance non-lineaire de 
zo avec 0 est compatible avec celle obtenus dans les simu¬ 
lations de particules frottantes IE10DI et a pour origine 
le desordre geometrique de ce type d’empilement Ena. 

La pression initiale Po augmente aussi au-dela du blo¬ 
cage a partir d’une pression residuelle sous 0j, qui comme 
pour zo, est due au figeage de la structure (Fig. ip 
Nous comprimons l’empilement de maniere non isotrope, 
en deplacant uniquement le piston mural. Malgre nos ef¬ 
forts pour preparer le materiau avec autant de precaution 
que faire se peut pn eng, en utilisant un protocole de 
compactification logarithmique sur une journee, l’empile- 
ment conserve de l’anisotropie clairement mis en evidence 
par l’existence d’une contrainte de cisaillement residuelle 
To proportionnelle a la pression Pq (Fig. HP- Cepen- 
dant, le rapport To/Po < 1, ce qui est attendu pour 
les empilements ou la contrainte de compression domine. 
Dans notre systeme, la sollicitation localisee est axi- 
symetrique. Par consequent, de telles fluctuations ani- 
sotropes m sont moyennees et notre mesure est robuste 
face a l’anisotropie. Il est a noter que la plupart des em¬ 
pilements prepares mimeriquement presentent aussi des 
contraintes de cisaillement residuelles et une forte aniso- 
tropie EH 02]. 


V. REPONSE AU GONFLEMENT 

Dans un premier temps, nous allons etudier le reseau 
de contact lors du gonflement de l’intrus. La figure [4] 
presente le changement de contact dans ce reseau pen¬ 
dant la phase de compression. On remarque que ce chan¬ 
gement de contact concerne plusieurs grains, repartis de 



Figure 3. Etat de contrainte initial (extraite de [32] 1. 

(a) Nombre moyen de contact initial zo (□); (b), pression 
Po (A) et contrainte de cisaillement To (v) en fonction de la 
fraction surfacique 0. La ligne noire est un a just ement a la 
fonction zo = z p ((f> — 0 j ) 0 ' 5 + zj, avec <f>j = 0.8251 ± 0.0009, 
z p = 10.0±0.5, et zj = 3.9=b0.1. La ligne en pointilles indique 
4>j- 

Initial stress state (adapted from \ 32$). (a) Initial average 
contact number zo ) ; (b), pressure Po (A) and shear stress 
To (V ) vs. packing fraction f>. The black line is a fit to zo = 
Zp^-cfjf-^Tzj, with <f)j = 0.8251P0.0009, = 10.0P0.5, 

and zj = 3.9 ±0.1. The dash line indicates <f>j . 
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maniere homogene dans l’empilement (fig. [4^i) . Ceci est 
d’autant plus remarquable que la solicitation elle-meme 
n’engendre pas une deformation homogene de l’empile- 
ment : le deplacement des grains est d’autant plus faible 
que Ton s’eloigne du centre de l’intrus (fig. [ 4 ) 3 ) . 

Ceci suggere fortement que la reponse mecanique est 
dominee par le desordre. Le nombre de changement de 
contact moyen par grain Az croit lineairement avec le 
parametre de dilatation a* (Fig. |4j^) , avec une pente po¬ 
sitive qui depend de <f>. Done, lors du gonflement de l’in- 
trus et du cisaillement de la structure, des contacts sont 
en moyenne crees, ce qui est en accord avec le fait que la 
surface tot ale occupee par les grains diminue, mais sans 
aucune structure spatiale ni correlation avec le forcage 




Figure 4. Augmentation du nombre de contacts. Car¬ 
tographies (a) des changements du nombre de contacts et (b) 
des deplacements radiaux u r pour <f> — 0.8294 et a* = 4.4%. 
Chaque grain colore en jaune a gagne ou perdu au moins un 
contact. Les grains sans couleurs sont masques par le tube 
pneumatique qui alimente l’intrus en air comprime. (c) : 
Difference du nombre de contact moyen avec l’etat initial, 
A z en fonction du parametre de gonflement a* = jjh pour 
differentes fractions surfaciques. Le code couleurs s’echelonne 
de bleu pour les faibles fractions surfaciques a rouge pour les 
grandes fractions surfaciques. (d) : Pente de la courbe A z- 
a* en fonction de la fraction surfacique, <f>. 

Increase of contact number. Maps (a) of contact changes 
and (b) of radial displacements u r for <f> — 0.8294 and 
a* — 4.4%. The yellow grains have lost or gained at least one 
contact. The uncolored grains sit below the pneumatic tube 
connected to the intruder, which masks the field of view, (c) 
Difference between the average contact number with the ini¬ 
tial state, A z vs. the inflation parameter a* = for various 
packing fractions. Color code spans from blue to red with in¬ 
creasing packing fractions, (d) Slope \z of the A z-a* curve 
vs. the packing fraction, (f. 


mecanique. Nous modelisons cet accroissement par une 
fonction lineaire et nous reportons la pente Xz dans la fi¬ 
gure W en fonction de <f>. Premierement, est toujours 
strictement positif. Cela suggere que de l’energie pour- 
rait etre dissipee lors d’un tel processus et que des effets 
non-lineaires sont attendus. Ensuite, la courbe presente 
un maximum precisement a </>j, qui est la fraction volu- 
mique de blocage. <f>j joue done un role precis dans la 
reponse mecanique de l’empilement. De telles observa¬ 
tions pourraient fournir de precieux indices pour la for¬ 
mulation d’un modele micro-mecanique. 

Dans la suite, nous allons considerer une formula¬ 
tion plus classique des tenseurs de contraintes et de 
deformations en utilisant leurs invariants. Nous allons 
considerer l’exces de pression P et de contrainte de ci¬ 
saillement r par rapport a l’etat de contrainte initial pour 
l’ensemble des a* et (f> mesures. Le comportement d’un 
materiau elastique lineaire, homogene et isotrope clas¬ 
sique a une relation lineaire entre l’etat de contrainte 
et l’etat de deformation, ce qui peut s’exprimer grace 
aux relations P = — Ke et r = 2Gy, ou K et G sont 
respectivement le module de compressibility et le mo¬ 
dule de cisaillement. Dans l’hypothese d’une geometrie 
2D axisymetrique, l’integration de l’equation d’equilibre 
conduit a r ~ G'j ~ a*/r 2 , ou r est la variable d’espace : 
la contrainte de cisaillement diminue avec la distance au 
centre de l’intrus. En revanche, la pression P ~ 1 ne 
depend pas de l’espace. Ainsi, dans cette formulation, 
le forcage localise induit des deformations volumiques et 
cisaillantes. 

Nous montrons en figure [5] les quatre cartes typiques 
des invariants pour une fraction d’empilement <f> = 0.8294 
superieure a <f>j et un taux de dilatation a* (4.4 x 10 -2 ). 
Nous observons de fortes fluctuations spatiales et la 
reponse s’ecarte de la reponse linaire elastique precisee 
plus tot. Cet ecart est inherent aux inhomogeneites 
presentes dans les materiaux desordonnes. Cependant, 
en considerant des echelles de longueurs plus grandes, on 
observe que l’axisymetrie du chargement est globalement 
respectee et que l’intensite de la reponse decroit avec la 
distance au centre de l’intrus. On peut ainsi faire l’hy- 
pothese que notre cellule est sufhsamment grande pour 
ne pas sentir les effets de confinement. Proche de l’in- 
trus, une import ante dilatation se produit a cause de 
la condition limite imposee par l’intrus lui-meme. Etant 
donne que l’intrus possede un diametre plus important 
que celui des grains, la variation de fraction surfacique 
localisee due au gonflement reste singuliere par rapport 
a la fraction surfacique moyenne. Le reste de l’empile- 
ment presente de large fluctuations (avec ici un ecart type 
de 3 x 10 -3 ) avec un leger accroissement de la fraction 
surfacique (d’amplitude 6 x 10- 5 ) (Fig.§0 , assurant la 
conservation globale du volume. 

Par consequent, en dehors de la premiere couronne de 
grains autour de l’intrus, que nous allons exclure de 1’ana¬ 
lyse, le materiau peut etre considere comme incompres¬ 
sible. Nous postulerons done dans la suite que e = 0. 
Ainsi, la solicitation mecanique se revele plutot etre 
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une deformation cisaillante pure. La deuxieme observa¬ 
tion importante montre un ecart significatif a la reponse 
elastique lineaire : le champ de pression est inhomogene et 
diminue a mesure que Ton s’ecarte du centre de l’intrus. 
Comme ce champ de pression ne correspond pas a des va¬ 
riations sensibles de volume, il est necessairement induit 
par la deformation de cisaillement. Cet effet est la mani¬ 
festation de la presence de dilatance dans un systeme a 
volume constant, un effet bien connu dans les milieux 
granulaires [29] . Le coefficient de dilatance a pression 
constante est defini par P = P 7 2 , et lie a celui defini 
a volume constant D par le module de compressibility K 
tel que R = DK [43] . 

Finalement, tandis que la pression P (moyennee ortho- 
radialement) evolue lineairement avec a* (Fig. Hi), 
la deformation de cisaillement 7 (moyennee ortho- 
radialement) croit plus vite avec a* (Fig. [BJd). C’est une 
indication de la nature non-lineaire du materiau. En re¬ 
vanche, le travail des efforts de cisaillement ry evolue en 
a * 2 (Fig. [Hjl) , ce qui suggere que la reponse mecanique 
du milieu granulaire a la sollicitation est dominee par des 
deformations elastiques. Ainsi, malgre le changement si¬ 
gnificatif de contacts dans l’empilement, il n’y a essentiel- 
lement pas de dissipation de l’energie. Ces observations 
sont valables pour l’ensemble des fractions surfaciques 


considerees. 


VI. RELATIONS CONSTITUTIVES 

Nous detaillons a present h analyse quantitative des re¬ 
lations de comport erne nt entre les invariants r, P et 7 . 
Nous rassemblons l’ensemble des donnees moyennees or- 
thoradialement P(r, a*) et r(r, a*) en fonction de 7 ( 7 *, a*), 
oil r est la distance au centre de l’intrus. Les figures [T^t 
et b montrent la contrainte de cisaillement et la pres¬ 
sion en fonction de la deformation de cisaillement 7 
pour differentes fractions surfaciques </>. Sous 0j, la pres¬ 
sion P et la contrainte de cisaillement r presented la 
dependance attendue avec la deformation de cisaillement 
7 : r = 2Go 7 et P — Po 7 2 , oil Go et Rq ne dependent 
pas de <fi. Au dessus de </>j, des non-linearites apparaissent 
sous la forme d’une diminution de la variation de r et P 
avec 7 . La transition de blocage joue un role essentiel et 
nous rutilisons pour decrire nos donnees avec un unique 
jeu de parametres definis par 

P =[R 0 + i? ra *(A</>, 7 )] 7 2 (10) 

t = 2 [G 0 + G n i(Acf>, 7 )] 7 ( 11 ) 



Figure 5. Cartographies des invariants des 
deformations et des contraintes (extraite de [32]). 
Cartographies de la dilatation, s ,(a), cisaillement, 7 , (b), 
pression, P, (c) et contrainte de cisaillement, r, (d), pour 
0 = 0.8294 et a* — 4.4 x 10 -2 . Les images representees ici ne 
represented qu’un tiers du systeme entier. 

Maps of the strain and stress invariants. (Adapted 
from [ Maps of dilation, e ,(a), shear strain, 7 , (b), 
pressure, P, (c) and shear stress, r, (d), for (f = 0.8294 and 
a* 4.4 x 10 -2 . The uncolored grains sit below the pneumatic 
tube connected to the intruder, which masks the field of view. 


Figure 6. Reponse au gonflement. Moyennes spatiales des 
invariants des tenseurs en fonction du parametre de gonfle¬ 
ment a*, (a) Pression P, (b) Cisaillement 7 , (c) Contrainte 
de cisaillement a et (d) puissance de travail des efforts de ci¬ 
saillement ycr. Le code couleurs s’echelonne de bleu pour les 
faibles fractions surfaciques a rouge pour les grandes fractions 
surfaciques. 

Response to inflation. Spatial averaged of the tensors vs. 
the inflation parameter a*. (a) Pressure P, (b) shear strain 7 , 
(c) shear stress a and (d) shear stress power ycr. Color code 
spans from blue to red with increasing packing fractions. 
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Figure 7. Relations de comportement (figure adaptee de 
32;). Pression, P (a), et contrainte de cisaillement, r (b), 
en fonction du cisaillement 7 , pour 21 fractions surfaciques 
4> G [0.8102 — 0.8343]. Les lignes pleines sont donnees par 
les equations . Le code couleurs s’echelonne de bleu 

pour les faibles fractions surfaciques a rouge pour les grandes 
fractions surfaciques. (c) et (d) : memes donnees que (a) et 
(b) normalisees par 7 C (0),P C (0) et r c (0). Les lignes pleines 
sont donnees par la version normalisee des equations, 
et les lignes en pointilles indiquent les regimes asymptotiques. 
Constitutive laws. (Adapted from f ~3^). Pressure, P (a), 
and shear stress, r (b), vs. shear strain, 7 , for 21 packing 
fractions G [0.8102 — 0.8343]. The solid lines are given by 
Eqs. fTdpTp . Color code spans from blue to red with increasing 
packing fractions, (c) and (d) : same data as (a) and (b) 
rescaled by 7 C (0),P C (0) and r c (<ft). The solid lines are given 
by the rescaled version of Eqs. and the dashed lines 

indicate the asymptotic regimes. 

avec A (f = (f - (fj, G 0 = 6.0 ± 0.2 x 1(T 2 , R 0 = 1.2 =b 
0.1 x 10 1 et 


Rni( A0,7) = 

G n /(A0,7) = 


0 for (j) < f)j 

aA0 M 7 a_2 for0>0j 

0 for <f> < <f>j 

bA(f> v ^~ x for <f> > <pj 


avec fi = 1.7d=0.1, a = 1.0A0.1, a = 8.1=b0.3x 10 2 , v = 
1.0 =b 0.1, p = 0.4 =b 0.1, b = 7.5 ± 0.3 x KT 1 . A partir 
de ces relations, on peut tracer deux courbes maitresses 
presentees dans les figures [ 7)3 et d avec 7 C ^ A0^, 
r c = 2Go^f c et P c = Ro7c- Bien que les exposants (/i, a) 
et (u, /?) aient ete obtenus independamment, il est remar- 
quable que ( = /i/(2 — a) et ( = v/{l — P) conduisent 
a la me me valeur ( = 1.7. Par consequent, les donnees 
de P et r presentent un comportement non-lineaire si- 
milaire et peuvent etre decrites de la meme maniere. Les 
equations et les lois d’echelles associees sont des 

resultats cles de l’etude. II est a noter que le regime 


lineaire observe ici ne doit pas etre confondu avec la 
reponse lineaire puisqu’elle se produit pour de grandes 
deformations. Nous expliquons la presence de ce nou¬ 
veau regime par une saturation des non-linearites. Nous 
pensons qu’un tel regime, qui existe aussi pour les em- 
pilements non bloques, merite une etude approfondie. 
Pour de petites deformations (7 cx 10 -6 ), sondees dans 
les etudes numeriques gl S3 mais plus faibles que les 
plus petites deformations sondees dans Lexperience, on 
s’attend a retrouver la veritable reponse lineaire pour 
tout A<f> > 0 [25]. Pour des deformations du meme 
ordre de grandeur que celles de Lexperience, des etudes 
numeriques recentes ont rapporte Lexistence d’un chan- 
gement de regime : pour des deformations plus grandes, 
la reponse devient non-lineaire, avec un exposant P ~ 
0.5 m sn , ce qui est compatible avec le regime non- 
lineaire observe ici. Dans de telles etudes, des particules 
avec un potentiel de deformation harmonique sont uti- 
lisees, i.e. Pq = A <f> et le croisement 7 * = Pq = A (j) b , ou 
b = 1. Ici, nos grains ne sont pas harmoniques et leur 
contact a un exposant b = 1.8. En outre, nous observons 
7 C ^ A<^, avec ( ~ b, ce qui induit 7 C 00 7 *. Ces deux 
seuils sont differents, mais ils ont le meme exposant, ce 
qui suggere qu’il se comportent de la meme maniere vis-a- 
vis de la transition de blocage. Cette etude preconise des 
investigations plus poussees, en particulier dans le regime 
non-lineaire, pour devoiler le role des non-linearites dans 
les lois de contact ainsi que celui de la friction entre 
grains. De plus, nous observons un role important de la 
dilatance dans nos grains cisailles, qui montrent que de 
tels effets non-lineaires doivent etre pris en compte pour 
decrire completement le materiau. L’importance de la di¬ 
latance dans les solides marginaux a recemment ete mise 
en lumiere dans [45^ , ou il a ete montre que le coefficient 
de Reynolds a volume constant Ry se comporte comme 
A0 -1 / 2 . Ici, nous reportons aussi un comportement sin- 
gulier, mais qui concerne les proprietes non-lineaires et 
en particulier badoucissement de la dilatance. Dans un 
contexte un peu different, Ren et al. m ont observe une 
tres forte augmentation de la dilatance sous cisaillement 
homogene lorsque la transition de blocage est approchee 
par le dessous. La valeur du coefficient de dilatance me- 
sure ici est tres grande (Rq ~ 10 4 N/m), et pourrait 
correspondre a une saturation de la divergence observee 
par Ren [31] . 


VII. PROFILS DE DEFORMATION 


Nous procedons a une verification de bauto-coherence 
du modele en injectant les relations constitutives du 
modele pour calculer les profils 7 et les comparer aux 
mesures experiment ales. La geometrie axisymetrique as¬ 
sure que la distance au centre r de bintrus est le seul 
parametre d’espace dont depend e et cr. Ainsi, en coor- 
donnees cylindriques, les relations constitutives peuvent 
s’ecrire sous la forme tensorielle 




















rj 

Figure 8 . Profils de cisaillement (extraite de [32] ). Profils 
de cisaillement pour (>) (0 = 0.8208; a* = 0.0374), (o) (0 = 
0.8268; a* = 0.0314) et (<) = 0.8338; a* = 0.0306). Les 

marqueurs correspondent aux donnees experimental es e t les 
lignes pleines a L integration numerique de l’equation |l2| . La 
ligne verte pointillee indique le seuil (r c ,7 c ) pour = 0.8268 
et a* = 0.0314. La zone grisee correspond a la zone occupee 
par le ” ballon”. 

Shear strain profiles (Adapted from J Ifflj). Shear strain pro¬ 
file for (>) ((/> = 0.8208; a* = 0.0374), (o) (</> = 0.8268; a* = 
0.0314; and (<) (f> — 0.8338; a* = 0.0306;. The symbols are 
experimental data and the solid lines come from the integra¬ 
tion of eg. \12l The green dashed line indicates the crossover 
for the case m 0.8268; a* = 0.0314;. 



Afin de simplifier, nous choisissons d’introduire la 
deformation de cisaillement reduite 7 = 7 / 7 c . La rela¬ 
tion precedente devient 


(T 


-p c {r+ 7) 


1 0 
0 1 


+ T c ( 7 ^ + 7 ) 


-1 0 
0 1 


L’equilibre mecanique se traduit par V-cr = 0, et conduit 
a 


--Pc («7 a 1 + 2 x) ~r~ Tc 77 1 + 1 )7 _2r< 


drf 


r+ 7 


dr 


dr 


= 0, 


qui, par separation des variables 7 et r, peut etre reecrit 
en 


P c (g7 q 1 + 27) + r c (/3 7^ 1 + 1) 

7^+7 


dS( ■ 


-2r c —. (12) 

r 


Nous integrons numeriquement l’equation [ 12 ] avec la 
condition limite 7 (r = rfi) = R*/ 7 C et nous obtenons 
les profils traces sur la figure [5] a laquelle nous avons 
ajoute les donnees experimentales. L’accord est satisfai- 
sant, etant donne qu’il n’y a aucun parametre ajustable 
et que nous avons neglige le confinement a grand r. Ceci 
confirme que notre description des relations constitutives 
est coherente. 


VIII. CONCLUSION 


Dans ce travail, nous avons fourni une caracterisation 
quantitative de la reponse elastique d’un empilement 2D 
de grains proche de la transition de blocage au gonfle- 
ment local d’un intrus. La geometric particuliere sonde la 
reponse a un cisaillement inhomogene a volume constant 
et peut etre vue comme un analogue a symetrie cylin- 
drique du probleme de l’inclusion d’Eshelby. Cette etude 
pourrait ainsi trouver une analogie avec les problemes de 
plasticite ou les contraintes elastiques sont redistributes 
par un evenement plastique localise mm- 

Nos resultats mettent en lumiere l’effet de dilatance 
et devoilent un regime non-lineaire au dela de la transi¬ 
tion de blocage oil le module de cisaillement et le mo¬ 
dule de dilatance diminuent jusqu’a ce qu’un nouveau 
regime lineaire apparaisse pour de grandes deformations. 
Cette etude montre l’existence d’un cisaillement critique 
7 C qui separe le regime non lineaire d’un regime lineaire 
sature. Son comportement avec la distance au blocage est 
coherent avec un autre croisement 7 *, recemment iden- 
tifie lorsque l’elasticite passe du regime lineaire au regime 
non-lineaire [43 , mais est tel que 7 * 7 c . Enfin, notre 
etude fournit les fondations pour un modele continu de 
tels materiaux. En particulier, l’observation du change- 
ment de reseau de contact est une observation cle pour 
pouvoir a terme obtenir une description micro-mecanique 
avec un formalisme de Cosserat [48^. 
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